Divergent streamlines and free vortices in Newtonian fluid flows in microfluidic flow focusing devices

The appearance of divergent streamlines and subsequent formation of free vortices in Newtonian fluid flows through microfluidic flow-focusing geometries is discussed in this work. The micro-geometries are shaped like a cross-slot but comprise three entrances and one exit. The divergent flow and subsequent symmetric vortical structures arising near the centreline of the main inlet channel are promoted even under creeping flow conditions, and are observed experimentally and predicted numerically above a critical value of the ratio of inlet velocities (VR). As VR is further increased these free vortices continue to grow until a maximum size is reached due to geometrical constraints. The numerical calculations are in good agreement with the experimental observations and we probe numerically the effects of the geometric parameters and of inertia on the flow patterns. In particular, we observe that the appearance of the central recirculations depends non-monotonically on the relative width of the entrance branches and we show that inertia enhances the appearance of the free vortices. On the contrary, the presence of the walls in three-dimensional geometries has a stabilizing effect for low Reynolds numbers, delaying the onset of these secondary flows to higher VR. The linearity of the governing equations for creeping flow of Newtonian fluids was invoked to determine the flow field for any VR as a linear combination of the results of three other independent solutions in the same geometry

The finite-volume method in computational rheology

The finite volume method (FVM) is widely used in traditional computational fluid dynamics (CFD), and many commercial CFD codes are based on this technique which is typically less demanding in computational resources than finite element methods (FEM). However, for historical reasons, a large number of Computational Rheology codes are based on FEM. There is no clear reason why the FVM should not be as successful as finite element based techniques in Computational Rheology and its applications, such as polymer processing or, more recently, microfluidic systems using complex fluids. This chapter describes the major advances on this topic since its inception in the early 1990’s, and is organized as follows. In the next section, a review of the major contributions to computational rheology using finite volume techniques is carried out, followed by a detailed explanation of the methodology developed by the authors. This section includes recent developments and methodologies related to the description of the viscoelastic constitutive equations used to alleviate the high-Weissenberg number problem, such as the log-conformation formulation and the recent kernel-conformation technique. At the end, results of numerical calculations are presented for the well-known benchmark flow in a 4:1 planar contraction to ascertain the quality of the predictions by this method

Steady flow of power-law fluids in a 1:3 planar sudden expansion

The laminar flow of inelastic non-Newtonian fluids, obeying the power-law model, through a planar sudden expansion with a 1:3 expansion ratio was investigated numerically using a finite volume method. A broad range of power-law indices in the range 0.2 n 4 was considered. Shear-thinning, Newtonian and shear-thickening fluids are analyzed, with particular emphasis on the flow patterns and bifurcation phenomenon occurring at high Reynolds number laminar flows. The effect of the generalized Reynolds numbers (based on power-law index, n, and the in flow channel height, h) on the main vortex characteristics and Couette correction are examined in detail in the range varying from 0.01 Regen 600. Values for the critical generalized Reynolds number for the onset of steady flow asymmetry and the appearance of a third main vortex are also included. We found that the shear-thinning behavior increases the critical Regen, while shear-thickening has the opposite effect. Comparison with available literature and with predictions using a commercial software (FluentR 6.3.26) are also presented and discussed. It was found that both results are in good agreement, but that our code is able to achieve converged solution for a broader range of flow conditions, providing new benchmark quality data

High performance microfluidic rectifiers for viscoelastic fluid flow

The flow of Newtonian and non-Newtonian fluids within microfluidic rectifiers with a hyperbolic shape was investigated to assess the effect of the bounding walls on the diodicity of the microfluidic device and achieve high flow anisotropy. Three microchannels were used, with different depths and the same geometrical configuration, which creates a strong extensional flow and generates high anisotropic flow resistance between the two flow directions. The Newtonian fluid, de-ionized water, was used as a reference fluid. The viscoelastic fluid used was an aqueous solution of polyethylene oxide (0.1% w/w) with high molecular weight. The flow patterns were visualized using streak photography and the velocity field was investigated using micro-particle image velocimetry. Moreover, pressure drop measurements were performed in order to compare the diodicity achieved in the microfluidic rectifiers. For the Newtonian fluid flow, the experimental results are compared with numerical predictions obtained using a finite-volume method and good agreement was found between both approaches. For the viscoelastic fluid, significant anisotropic flow resistance can be achieved. The effect of the bounding walls was analysed and found to be qualitatively similar for all microchannels. Nevertheless, in quantitative terms, the diodicity is enhanced when the wall effect is reduced, i.e. when the channels are deeper. A maximum diodicity above six was found for the deeper channel, a value well beyond those previously reported

Laminar flow in three-dimensional square-square expansions

In this work we investigate the three-dimensional laminar flow of Newtonian and viscoelastic fluids through square–square expansions. The experimental results obtained in this simple geometry provide useful data for benchmarking purposes in complex three-dimensional flows. Visualizations of the flow patterns were performed using streak photography, the velocity field of the flow was measured in detail using particle image velocimetry and additionally, pressure drop measurements were carried out. The Newtonian fluid flow was investigated for the expansion ratios of 1:2.4, 1:4 and 1:8 and the experimental results were compared with numerical predictions. For all expansion ratios studied, a corner vortex is observed downstream of the expansion and an increase of the flow inertia leads to an enhancement of the vortex size. Good agreement is found between experimental and numerical results. The flow of the two non-Newtonian fluids was investigated experimentally for expansion ratios of 1:2.4, 1:4, 1:8 and 1:12, and compared with numerical simulations using the Oldroyd-B, FENE-MCR and sPTT constitutive equations. For both the Boger and shear-thinning viscoelastic fluids, a corner vortex appears downstream of the expansion, which decreases in size and strength when the elasticity of the flow is increased. For all fluids and expansion ratios studied, the recirculations that are formed downstream of the square–square expansion exhibit a three-dimensional structure evidenced by a helical flow, which is also predicted in the numerical simulations

Effect of the contraction ratio upon viscoelastic fluid flow in three-dimensional square-square contractions

In this work we investigate the laminar flow through square–square sudden contractions with various contraction ratios (CR¼2.4, 4,8and12), using a Newtonian fluid and a shear-thinning viscoelastic fluid. Visualizations of the flow patterns were carried out using streakline photography and detailed velocity field measurements were performed using particle image velocimetry. The experimental results are compared with numerical predictions obtained using a finite-volume method. For the Newtonian fluid, a corner vortex is found upstream of the contraction and increasing flow inertia leads to a reduction of the vortex size. Good agreement is observed between experiments and numerical simulations. For the shear-thinning fluid flow a corner vortex is also observed upstream of the contraction independently of the contraction ratio. Increasing the elasticity of the flow, while still maintaining low inertia flow conditions, leads to a strong increase of the vortex size, until an elastic instability sets in and the flow becomes time-dependent at DeE200, 300, 70 and 450 for CR¼2.4, 4, 8 and 12, respectively. At low contraction ratios, viscoelasticity brings out an anomalous divergent flow upstream of the contraction. For both fluids studied the flow presents a complex three-dimensional helical vortex structure which is well predicted by numerical simulations. However, for the viscoelastic fluid flow the maximum Deborah number achieved in the numerical simulations is about one order of magnitude lower than the critical Deborah number for the onset of the elastic instability found in the experiments

Three-dimensional flow of Newtonian and Boger fluids in square-square contractions

The flow of a Newtonian fluid and a Boger fluid through sudden square–square contractions was investigated experimentally aiming to characterize the flow and provide quantitative data for benchmarking in a complex three-dimensional flow. Visualizations of the flow patterns were undertaken using streakline photography, detailed velocity field measurementswere conducted using particle image velocimetry (PIV) and pressure drop measurements were performed in various geometries with different contraction ratios. For the Newtonian fluid, the experimental results are compared with numerical simulations performed using a finite volume method, and excellent agreement is found for the range of Reynolds number tested (Re2 ≤23). For the viscoelastic case, recirculations are still present upstream of the contraction but we also observe other complex flow patterns that are dependent on contraction ratio (CR) and Deborah number (De2) for the range of conditions studied: CR = 2.4, 4, 8, 12 and De2 ≤150. For low contraction ratios strong divergent flow is observed upstream of the contraction, whereas for high contraction ratios there is no upstream divergent flow, except in the vicinity of the re-entrant corner where a localized a typical divergent flow is observed. For all contraction ratios studied, at sufficiently high Deborah numbers, strong elastic vortex enhancement upstream of the contraction is observed, which leads to the onset of a periodic complex flow at higher flow rates. The vortices observed under steady flow are not closed, and fluid elasticity was found to modify the flow direction within the recirculations as compared to that found for Newtonian fluids. The entry pressure drop, quantified using a Couette correction, was found to increase with the Deborah number for the higher contraction ratios

Flow of a blood analogue solution through microfabricated hyperbolic contractions

The flow of a blood analogue solution past a microfabricated hyperbolic contraction followed by an abrupt expansion was investigated experimentally. The shape of the contraction was designed in order to impose a nearly constant strain rate to the fluid along the centerline of the microgeometry. The flow patterns of the blood analogue solution and of a Newtonian reference fluid (deionized water), captured using streak line imaging, are quite distinct and illustrate the complex behavior of the blood analogue solution flowing through the microgeometry. The flow of the blood analogue solution shows elastic-driven effects with vortical structures emerging upstream of the contraction, which are absent in Newtonian fluid flow. In both cases the flow also develops instabilities downstream of the expansion but these are inertia driven. Therefore, for the blood analogue solution at high flow rates the competing effects of inertia and elasticity lead to complex flow patterns and unstable flow develops

Extensional rheology and elastic instabilities of a wormlike micellar solution in a microfluidic cross-slot device

Wormlike micellar surfactant solutions are encountered in a wide variety of important applications, including enhanced oil recovery and ink-jet printing, in which the fluids are subjected to high extensional strain rates. In this contribution we present an experimental investigation of the flow of a model wormlike micellar solution (cetyl pyridinium chloride and sodium salicylate in deionised water) in a well-defined stagnation point extensional flow field generated within a microfluidic cross-slot device. We use micro-particle image velocimetry (m-PIV) and full-field birefringence microscopy coupled with macroscopic measurements of the bulk pressure drop to make a quantitative characterization of the fluid’s rheological response over a wide range of deformation rates. The flow field in the micromachined cross-slot is first characterized for viscous flow of a Newtonian fluid, and m-PIV measurements show the flow field remains symmetric and stable up to moderately high Reynolds number, Re z 20, and nominal strain rate, _3nom z 635 s1. By contrast, in the viscoelastic micellar solution the flow field remains symmetric only for low values of the strain rate such that _3nom # lM1, where lM ¼ 2.5 s is the Maxwell relaxation time of the fluid. In this stable flow regime the fluid displays a localized and elongated birefringent strand extending along the outflow streamline from the stagnation point, and estimates of the apparent extensional viscosity can be obtained using the stressoptical rule and from the total pressure drop measured across the cross-slot channel. For moderate deformation rates (_3nom $ lM1) the flow remains steady, but becomes increasingly asymmetric with increasing flow rate, eventually achieving a steady state of complete anti-symmetry characterized by a dividing streamline and birefringent strand connecting diagonally opposite corners of the cross-slot. Eventually, as the nominal imposed deformation rate is increased further, the asymmetric divided flow becomes time dependent. These purely elastic instabilities are reminiscent of those observed in crossslot flows of polymer solutions, but seem to be strongly influenced by the effects of shear localization of the micellar fluid within the microchannels and around the re-entrant corners of the cross-slot

Flow of low viscosity Boger fluids through a microfluidic hyperbolic contraction

In this work we focus on the development of low viscosity Boger fluids and assess their elasticity analyzing the flow through a microfluidic hyperbolic contraction. Rheological tests in shear and extensional flows were carried out in order to evaluate the effect of the addition of a salt (NaCl) to dilute aqueous solutions of polyacrylamide at 400, 250, 125 and 50 ppm (w/w). The rheological data showed that when 1% (w/w) of NaCl was added, a significant decrease of the shear viscosity curve was observed, and a nearly constant shear viscosity was found for a wide range of shear rates, indicating Boger fluid behavior. The relaxation times, measured using a capillary break-up extensional rheometer (CaBER), decreased for lower polymer concentrations, and with the addition of NaCl. Visualizations of these Boger fluids flowing through a planar microfluidic geometry containing a hyperbolic contraction, which promotes a nearly uniform extension rate at the centerline of the geometry, was important to corroborate their degree of elasticity. Additionally, the quantification of the vortex growth upstream of the hyperbolic contraction was used with good accuracy and reproducibility to assess the relaxation time for the less concentrated Boger fluids, for which CaBER measurements are difficult to perform